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Abstract 

Leading-log results are derived for the shear viscosity, electrical conductivity, 
and flavor diffusion constants in both Abelian and non-Abelian high temper- 
ature gauge theories with various matter field content. 

I. INTRODUCTION AND SUMMARY 

Transport coefficients, such as viscosities, diffusivities, or conductivity, characterize the 
dynamics of long wavelength, low frequency fluctuations in a medium. In condensed matter 
applications transport coefficients are typically measured, not calculated from first princi- 
ples, due to the complexity of the underlying microscopic dynamics. But in a weakly coupled 
quantum field theory, transport coefficients should, in principle, be calculable purely the- 
oretically. Knowledge of various transport coefficients in high temperature gauge theories 
is important in cosmological applications such as electroweak baryogenesis as well as 
hydrodynamic models of heavy ion collisions ||. 

In this paper, we consider the evaluation of transport coefficients in weakly coupled high 
temperature gauge theories. "High temperature" is taken to mean that the temperature 
is much larger than the zero-temperature masses of elementary particles, and any chemical 
potentials. In QED, this means T 3> m e , while in QCD, we require both T ^> Aqcd and 
T 3> m q . Corrections suppressed by powers of temperature (m q /T, Aq CD /T, etc.) will be 
ignored. This means that each transport coefficient will equal some power of temperature, 
trivially determined by dimensional analysis, multiplied by some function of the dimension- 
less coupling constants of the theory. 

As an example, the shear viscosity in a (single component, real) |y0 4 scalar theory has 
the high temperature form 



with a = 3033.5, up to relative corrections suppressed by higher powers of A The 
leading 1/A 2 behavior reflects the fact that the two-body scattering cross section is 0(A 2 ), 
and that transport coefficients are inversely proportional to scattering rates. 

In gauge theories, the presence of Coulomb scattering over a parametrically wide range of 
momentum transfers (or scattering angles) causes transport coefficients to have a more com- 
plicated dependence on the interaction strength. In QED, for example, the high temperature 
shear viscosity has the form 

T 3 

V = k—. r , (1.2) 

cr ma 1 

up to relative corrections which are suppressed by additional factors of 1/ In a -1 . Evaluating 
the overall constant k, while ignoring all terms suppressed by additional powers of 1/ lna -1 
(or powers of a) amounts to a "leading-log" calculation of the transport coefficient. 

In this paper, we will present leading-log calculations of the shear viscosity, electrical 
conductivity, and flavor diffusion constants in high temperature gauge theories (Abelian or 
non-Abelian) with various matter field content .0 It should be emphasized, however, that 
these leading-log results cannot be presumed to provide a quantitatively reliable determi- 
nation of transport coefficients in any real application. Gauge couplings in the standard 
model are never so tiny that corrections suppressed by 1/ In a" 1 , 1/ In a" 1 , or even 1/ lna^ 
are negligibly small. Nevertheless, the leading-log analysis of transport coefficients is a use- 
ful first step. In a companion paper, we extend our treatment and obtain "all-log" results 
which include all terms suppressed only by inverse logarithms of the gauge coupling (but 
drop sub-leading effects suppressed by powers of the coupling) [|]. 

Previous efforts to determine transport coefficients in hot gauge theories include many 
applications of relaxation time approximations [|7|-[T6[|, in which the full momentum depen- 



dence of relevant scattering rates are crudely characterized by a single relaxation time. Such 
treatments can, at best, obtain the correct leading parametric dependence on the coupling 



and a rough estimate of the overall coefficient (though some of them [|I3| -|T6| do not obtain 
the right parametric behavior). In addition, there have been a number of papers reporting 
genuine leading- log evaluations of various transport coefficients JTT|-|2"3"|] . However, we find 
that almost all of these results are incorrect due to a variety of both conceptual and technical 
errors. (In some cases ||17|| , the errors are numerically quite small.) For each transport coef- 
ficient we consider, specific comparisons with previous work will be detailed in the relevant 
section below. 

In section |T| we discuss how one may construct a linearized kinetic theory which is 
adequate for computing correctly the transport coefficients we consider, up to corrections 
suppressed by powers of coupling. We also show how the actual calculation of a transport 



1 The value quoted for a actually comes from our own evaluation of the (j) 4 shear viscosity, using the 
variational formulation described in section ||. This allows a higher precision evaluation than that obtained 
by discretizing the requisite integral equation as described in pi . 

2 We will not consider the bulk viscosity. It requires a significantly different, and more complicated, analysis 
than other transport coefficients. We also do not treat thermal conductivity, which is not an independent 
transport coefficient in the absence of nonzero conserved charges (besides energy and momentum). 
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coefficient may be converted into a variational problem; this is very convenient for numerical 
purposes. (Related, but somewhat different variational formulations appear in the litera- 
ture.) For leading- log calculations one may greatly simplify the resulting collision integrals, 
since the coefficient of the leading-log is only sensitive to small angle scattering processes. 
This is discussed in section PTT| . The details of the analysis for the electrical conductivity, 
flavor diffusivities, and shear viscosity are presented in sections |IV|, |V|, and |Vj, respectively. 
Throughout this paper, we will present results for arbitrary simple gauge groups, rather than 
specializing to SU(iV). Our notation for group factors (Cr, T r , g? r ), their SU(iV) values, 
and how to interpret them for Abelian problems, is explained in Appendix [A]. 

In the remainder of this introduction, we review the basic definitions of the various 
transport coefficients, and then summarize our results. 

A. Definitions 

At sufficiently high temperature, the equilibrium state of any relativistic field theory may 
be regarded as a relativistic fluid. The stress energy tensor T^ v defines four locally-conserved 
currents whose corresponding conserved charges are, of course, the total energy and spatial 
momentum of the system. At any point in the system, the local fluid rest frame is defined 
as the frame in which the local momentum density vanishes^] 

T°j(a;) = -v=> [local fluid rest frame at x] . (1.3) 

If the fluid is slightly disturbed from equilibrium, then the non-equilibrium expectation 
of Tjjj,, in the local fluid rest frame, satisfies the constitutive relation^ 

(Tij) = Sij (^)-T) [V< Uj + Vj m-l S id V' Ul ] - C 6^ V' u\ , (1.4) 

together with the exact conservation law <9 M (T Ml/ ) = 0. In the constitutive relation ( |1.4|) , rj is 
the shear viscosity, ( is the bulk viscosity, and u is local flow velocity. For small departures 
from equilibrium, 

u i = {T° i )/{E+^), (1.5) 

where e = T 00 is the energy density and s = \T l i the (local) pressure. The combination 
e + & is also known as the enthalpy. The constitutive relation ( |1.4| ) holds up to corrections 
involving further gradients or higher powers of u. 

In a similar fashion, in any theory containing electromagnetism (i.e., a U(l) gauge field), 
the electric current density j™ is conserved (d^j™ = 0) and satisfies the constitutive 
relation 



3 This is often termed the Landau-Lifshitz convention. In theories with additional conserved currents, such 
as a baryon number current, one may alternatively define the local rest frame as the frame in which there 
is no baryon number flux. This is the Eckart convention. As we wish to consider, among others, pure gauge 
theories in which no other conserved currents are present, the Landau-Lifshitz convention provides the only 
uniform definition of local flow. 

4 We use (—+++) metric conventions. 
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(1.6) 



for small departures from equilibrium. Here a is the (DC) electrical conductivity, 

And in theories containing one or more conserved "flavor" currents, j", which are not 
coupled to dynamical gauge fields (such as baryon number or isospin currents in QCD), 
these currents will satisfy diffusive constitutive relations, 

O'f ) = ~D aP Viin') , (1.7) 

where n a = (j a )° is a conserved charge density, and D = \ \D at3 \ \ is, in general, a matrix of 
diffusion constants. Once again, the constitutive relation (|1.7|) holds in the local fluid rest 
frame; otherwise, an additional convective (n a ) Vi term is present. 

In the limit of arbitrarily small gradients (so that the scale of variation in (T^) or (j^) 
is huge compared to microscopic length scales) and arbitrarily small departures from equi- 
librium, the constitutive relations (|1.4[), (|1.6j), and ( |1.7| ) may be regarded as definitions of 
the shear and bulk viscosities, electrical conductivity, and flavor diffusion constants. Alter- 
natively, one may use linear response theory to relate non-equilibrium expectation values 
to equilibrium correlation functions ||. This leads to well-known Kubo relations which 
express transport coefficients in terms of the zero-frequency slope of spectral densities of 
current-current, or stress tensor-stress tensor correlation functions, 



V = — lim l - / d 4 x e wt ([7r, m (t, x), vr im (0)]) eq , (1.8a) 

C= \ lim - / d A xe^ <[£>(*, x), ^(0)]) eq , (1.8b) 

a= l^I ^? U &X-)> jf M (0)]) eq , (1.8c) 

D aP = l^U d 4 x e*"* (\j?(t, x), j'7(0)]) eq . (1.8d) 



In Eq. ( |1.8a| ), 7r i?n = T i?n — Si m ^ denotes the traceless part of the stress tensor, while in 
Eq. (|1.8d|) , H = ||H a/ g|| is the "charge susceptibility" matrix describing mean-square global 



charge fluctuations (per unit volume), 

(1.9) 



■ a ^^ = l[(N«Ne)-(N«)(NP) 



where N a = j (<i 3 x) n a is a conserved charge, is the spatial volume, and (3 is the inverse 
temperature. 



B. Results 

1. Electrical conductivity 
The high temperature electrical conductivity has the leading-log form 

o = C T (1.10) 

e 2 me 1 
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TABLE I. Leading-log conductivity for various numbers of lcptonic charge carriers (N\ cptons ) and effec- 
tive number of leptonic plus quark scatterers (N spoc i os ) . The number A^ spec i es is a sum over all fermion fields 
weighted by the square of their electric charge. 



where the dimensionless coefficient C depends on the number (and relative charges) of the 
electrically charged matter fields which couple to the photon. The T/(e 2 hie -1 ) dependence 
of cr can be qualitatively understood as arising from the form a ~ e 2 T 2 r, where r ~ 
[e 4 T loge -1 ] 1 is the characteristic time scale for large-angle scattering (from either a single 
hard scattering or a sequence of small-angle scatterings which add up to produce a large 
angle deflection), or "transport mean free time."(] This is just the classic Drude model form 
a ~ ne 2 T/m, appropriately generalized to an ultra-relativistic setting (so n ~ T 3 , and m 
is replaced by the typical energy T). Our exact results for the leading- log coefficient C, for 
various subsets of the fermions of the standard model, are shown in table |. Which entry in 
the table to use depends on the temperature; each entry is valid when the included fermions 
are light (m -C T) but the excluded fermions are heavy. 

The dependence on matter field content is not precisely given by any simple analytic 
formula. But as detailed in section 0, the conductivity is approximately equal to 

j ^/ l2 4 C(3) 2 vr-3iV lcptons \ T 

V 37r 2 + 32iV species J e 2 \ne~ 1 ' V ' ; 

where A^i eptons is the number of leptonic charge carriers (not counting anti-particles separately 
from particles), and iV S p ec i es is a sum over all (Dirac) fermion fields weighted by the square 
of their electric charge assignments. So each leptonic species contributes 1 to N spec i eB , each 
down type quark contributes 1/3 [which is (— 1/3) 2 times 3 colors], and each up type quark 
contributes 4/3. It may be noted that this is exactly the same sum over charged species 
which appears in the (lowest-order) expression for the high temperature Debye screening 
mass for the photon, 

m 2 D = \ N species e 2 T 2 . (1.12) 



The approximate form ( |1.11|) is the result of a one-term variational approximation which 



becomes exact if the correct quantum statistics is replaced by classical Boltzmann statistics. 



5 The transport mean free time is the inverse collision rate for large-angle scattering (i.e., an O(l) change 
in direction). In high temperature gauge theories, the mean free time for any scattering is dominated by 
very small angle scattering and is of order (g 2 T) -1 , up to logarithms, while the transport mean free time is 
order (g 4 T 'In g^ 1 )' 1 @,| (jjjijllfff] . 
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This approximation reproduces the true leading-log coefficient to within an accuracy of 
better than 0.4% for all cases shown in Table |. In practice, neglect of subleading effects 
suppressed by powers of 1/ In e -1 is a far larger issue than the accuracy with which Eq. (|1.11 ) 
reproduces the exact leading- log coefficient. 

Our Eq. ( |1 . 1 1| ) is similar to the expression found by Baym and Heiselberg [^UJ. The most 
significant difference is that their expression is missing the 37r 2 term in the denominator. 
This term arises from Compton scattering and annihilation to photons — processes neglected 
in Ref. 0. 

At temperatures above the QCD scale, where scattering from quarks must be included, 
these leading-log results neglect the contribution of quarks to the electric current density j™ 
itself. This is quite a good approximation since the rate of strong interactions among quarks 
is much greater than their electromagnetic interactions, and will wash out departures from 
equilibrium in quark distributions much faster than the relaxation of fluctuations in lepton 
distributions, which depends on electromagnetic interactions. This simplification amounts 
to the neglect of corrections to the coefficient C suppressed by a| M /a;?. There are also 
relative 0(a s ) corrections arising from QCD effects when a lepton scatters from a quark. 

Weak interactions are ignored altogether in these calculations, so the results above, when 
applied to the standard model, are relevant at temperatures small compared to M w , but 
large compared to the masses of the quarks and leptons considered. In particular, the 
electromagnetic mean free path for large angle scattering must be small compared to the 
weak interaction mean free path which, parametrically, requires that a 2 3> a^(T/M^) 4 . 

For temperatures comparable or large compared to Mw, where the electroweak sector of 
the standard model is in its "unbroken" high temperature phase,(] the dynamics of the U(l) 
hypercharge field may be characterized by a hypercharge conductivity in complete analogy 
to ordinary electromagnetism. Neglecting relative corrections suppressed by tan 4 #\y and 
a'/a s , and negligible effects due to Yukawa couplings of right-handed leptons, the hyper- 
charge conductivity is determined by the hypercharge-mediated scattering of right-handed 
leptons. (Quarks and left-handed leptons scatter much more rapidly due to SU (2) or SU(3) 
interactions, and hence contribute much less to the conductivity than right-handed leptons.) 
The appropriate generalization of Eq. ( |1 . 1 1| ) has exactly the same form but with the charge e 
replaced by the hypercharge coupling g', N\ cptons replaced by half the number of right-handed 
leptons, and iV S p ec i es now given by the sum of the square of the hypercharge of each complex 
scalar field, plus half the square of the hypercharge of each chiral fermion field. For the 
three generations of the standard model, this means iV le p toIls = 3/2, and N species = 5 + n s /2, 
where n s is the number of Higgs doublets,[] so that leading-log hypercharge conductivity is 



6 Depending on details of the scalar sector, there need not be any sharp electroweak phase transition. The 
"high temperature electroweak phase" should be understood as the regime where the effective mass of the 
weak gauge bosons comes predominantly from thermal fluctuations, not the Higgs condensate. In other 
words, T <; v(T), where v(T) is the (temperature dependent) Higgs expectation value. 

7 We normalize "h yper charge" Y by Q = T3 + Y (as opposed to the other convention that Q = T3 + Y/2). 
The conductivity (1.13|) does not depend on this normalization convention. In greater detail, N Bpcc [ cs — 



2n s (l/2) 2 + 3/2 [(-l) 2 + 2(-l/2) 2 + 3(2/3) 2 + 3(-l/3) 2 + 6(l/6) 2 ] , where the various terms in the bracket 
come from right-handed leptons, left-handed leptons, right-handed up type quarks, right-handed down type 
quarks, and left-handed quarks, respectively. The Debye mass of the U(l) hypercharge gauge field again 



6 



Nf rj x (g 4 /T 3 ) lng" 1 

~0 27126 

1 60.808 

2 86.473 

3 106.664 

4 122.958 

5 136.380 

6 147.627 



TABLE II. Leading-log shear viscosity as a function of the number of (fundamental representation) 
fermion flavors with m<T, for gauge group SU(3). 



(approximately) 



a hyper ~6 4 C(3) 2 tt- 3 



?H 20 2 

1 1 — n s 

8 3 3 



T 



g/2 J n gl 1 



;i.i3) 



2. Shear viscosity 

The high temperature shear viscosity in a gauge theory with a simple gauge group (either 
Abelian or non-Abelian) has the leading-log form 



T 3 



7] = K 



g 4 in^- 1 



;i.i4) 



where g is the gauge coupling. For the case of SU(3) gauge theory (i.e., QCD), our results 
for the leading-log shear viscosity coefficient k for various numbers of fermion species are 
shown in table |TJ. 

The analysis may be easily generalized to an arbitrary gauge group with Nf Dirac 
fermions in any given representation. Once again, the numerical results are approximately 
reproduced by a relatively simple analytic form which is the result of a one term variational 
calculation, 



77~ 270d A C(5r - (v ] c- Y v) 



71 



T 3 



g 4 in^- 1 



;i.i5) 



where c is the 2x2 matrix 



c = (d A C A + Nf d F C F ) 



d A C A 
i 



jNf d F C F 



9tt^ 
128 



N t d F Cp d A 



1 -1 
1 1 



(1.16) 



and 



satisfies Eq. (1.12), with 



rj2 



d A 

f «F 



4N f d 



;i.i7) 



Here, dp and Cp denote the dimension and quadratic Casimir of the fermion representation, 
while g?a and C\ are the dimension and Casimir of the adjoint representation. (See Appendix 
[A].) In all cases studied, the expression ( |1.15[ ) is accurate to within 0.7%. The two-by-two 



matrix structure of expression ( |1.15| ) arises from the fact that the leading-log shear viscosity 
is sensitive to all two-particle scattering processes: fermion-fermion, fermion-gluon, and 
gluon-gluon. In particular, the non-diagonal second term in Eq. ( |1.16|) arises from Compton 



scattering and qq annihilation to gluons, as will be described in sections |T| and |VT 



An earlier result of Baym, Monien, Pethick, and Ravenhall JT7| coincides with 



Eqs. (|1.15|) - (|1.17|) except that they omitted this second term in the matrix c. A later 



paper of Heiselberg jlq] essentially agrees with our leading-log result for the shear viscosity 
of pure gauge theory,^ but the later treatment of fermions in this paper also missed the 
Compton scattering and annihilation contributions, and made additional errors not present 

0- 



in 



Plugging C*a = and d\ = dp = Cf = 1 into Eq. ( |1.15| ) yields (a good approximation to) 
the leading-log QED result for Nf charged leptons and no quarks, which is again accurate 
to within 0.7%. For an e + e~ plasma (iVf = 1), this gives 

w n9 /2\ 5 /529 128\ T 3 T 3 . 

7?~ 270C5 2 - + 7 = 1 87.129— -. 1.18 

' KJ \ttJ V112 oWeMne" 1 e 4 lne" 1 v ' 

The complete leading-log calculation for this case gives rj = 188.38 T 3 /(e 4 lne" 1 ). 

When applied to the standard model at temperatures small compared to Mw (but large 
compared to Aq CD ), the results of Table [TJ or Eq. ( |1.15| ), with g the running QCD coupling 
(evaluated at a scale of order T), are relevant for hydrodynamic fluctuations in a quark- 
gluon plasma occurring on length scales which are large compared to the strong interaction 
transport mean free path of order (g^T In g^ 1 )^ 1 , but small compared to the electromagnetic 
transport mean free path of order (e 4 T In e -1 ) -1 . In this regime, leptons may be regarded 
as freely streaming and decoupled from the quark-gluon plasma. 

On longer length scales, large compared to the electromagnetic transport mean free path, 
the shear viscosity is dominated by electromagnetic scatterings of out-of-equilibrium charged 
leptons. [Photons do not contribute significantly because they are thermalized by — > qq 
and 7g — > qg processes, whose rates are 0(aEM«s) and hence rapid compared to purely 
electromagnetic scatterings.] In this domain, the shear viscosity is approximately given by 



This form reproduces the correct leading-log coefficient to within 0.5%, but neglects strong 
interaction effects which give relative corrections suppressed by a s or oem/c^s (as well as 



8 The sign of the difference between the one-term ansatz result and the correct leading-log coefficient is 
reported incorrectly in Appendix A of Ref. In a related matter, Ref. (l8) incorrectly asserts that the 

exact 77 is a minimum of the variational problem set up in that paper; it is actually a maximum. 
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next-to-leading log corrections formally down by l/lna^). Once again, A^ eptons is the 
number of light leptonic species, and A^ S p ecies is the sum over all light fermion fields weighted 
by the square of their electric charges. 

The result ( |1.19| ) assumes that neutrinos may still be regarded as freely streaming and 
decoupled, so it is valid only on length scales small compared to the neutrino mean free 
path, which is of order [M z /T) 4 / '(aj^T). On scales large compared to the neutrino mean 
free path, the shear viscosity is dominated by neutrino transport and scales as T 4 times this 
mean free path, 

n = o(m. (1.20) 

Ww T ) 

We have not calculated the precise coefficient and, to our knowledge, no quantitative calcu- 
lation of neutrino viscosity is available in the literature.^ 

Finally, in the high temperature electroweak phase the viscosity, like the hypercharge 
conductivity, is dominated by right handed lepton transport. Neglecting relative order 
tan 4 W and a! /a s corrections plus negligible Yukawa coupling effects, the leading-log shear 
viscosity in this regime is given by Eq. ( |1.19| ) with e replaced by g', iVi eptons — > 3/2, and 

-^species * 5 -|- ?7. s /2. 

3. Baryon and lepton number diffusion 

To leading-log order, the fermion number diffusion constant in a QED or QCD-like gauge 
theory has the form 

D F = A F J — , (1.21) 

where g is the gauge coupling and A F is a constant. For the case of SU(3) gauge theory, 
our results for the fermion number diffusion coefficient A F for various numbers of fermion 
flavors are shown in table |TJ. 

As will be discussed in detail below, the diffusion constant depends on the rate at which 
a fermion scatters off either another fermion, or a gluon (photon) present in the high tem- 



perature plasma. In Table [TJ, the line labeled "0" is analogous to a quenched (or valence) 
approximation, and shows the result when only scattering off thermal gluons is included. 
The fairly weak dependence of the coefficient on Nf shows that in all cases the gluonic 
scattering contribution is dominant. 

As in the previous cases, these results for an SU(3) gauge theory are approximately 
reproduced by the simple analytic form 

/2 4 3 6 C(3) 2 tt- 3 \ T~ l 



24 + 4iV f +7T 2 J g A \ng- 1 



(1.22) 



3 A fairly careful estimate of the neutrino mean free path has been given in Ref. |27| 
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Nf D F x(g 4 T)\ng- 

"0" 16.0597 

1 14.3677 

2 12.9990 

3 11.8688 

4 10.9197 

5 10.1113 

6 9.4145 



TABLE III. Leading-log diffusion constant for fermion number density as a function of the number 
of fermion flavors with m < T in SU (3) gauge theory. The line labeled "0" shows the result when only 
scattering off thermal gluons is included. 



which is the result of a one-term variational approximation. In all cases studied this ex- 
pression is accurate to within 0.3%. The three factors in the denominator arise, in order, 
from t-channel gluon exchange with a gluon, i-channel gluon exchange with a quark, and 
Compton scattering or annihilation to gluons. 

This expression can be generalized to arbitrary simple gauge group and matter fields in 
any representation. The leading-log diffusion constant for the net number density of fermion 
flavor a is (approximately) given by 



n ~6«C(3) S 

IS n 



fa 3tt 2 

Y2 T Rb^b + — c Ra 



b 



T -l \ 

(1.23) 



where the sum is over all flavors and helicities b of the excitations that fermion a can scatter 
from in ab —>■ ab processes mediated by gauge-boson exchange, including separately particles 
and antiparticles. (So 2 terms appear for scattering off of a gauge boson, complex scalar, 
or Weyl fermion, and 4 terms for scattering from a Dirac fermion.) We have introduced the 
notation "ffh" over the sum as a reminder that the sum includes flavors [f], anti-flavors if 
distinct [f], and helicities [h]. The group representation normalization factor Tr = Crc^r/c^a 
is defined in appendix |A]. The symbol A& is 

. f 2, if b is a boson; ,^ 
b \ 1, if b is a fermion. 

If the relevant scattering is by photon exchange, then d& = d,R = 1, and e 2 a = g 2 C Ra is the 
squared electric charge of species a. The result ( |1 . 23| ) neglects any Yukawa interactions with 
scalar fields. Once again, the sum appearing in Eq. (|1.23p also appears in the lowest-order 
expression for the high temperature Debye mass, now generalized to an arbitrary simple 
gauge group and arbitrary matter content, 




m D = 2_J T R \ b g T . (1.25) 



When applied to the standard model at temperatures small compared to Mw, the result 
( |1.22j ) with g 2 = g 2 gives (a good approximation to) the leading-log result for the diffusion 
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constant which is appropriate for describing relaxation of fluctuations in baryon density, or 
equivalently the net density of any particular quark flavor, on scales which are large compared 
to the strong interaction transport mean free path of order (g 4 Tln g^ 1 )^ 1 . If the departure 
from equilibrium of the various quark densities has vanishing electric charge density, then the 
resulting relaxation is purely diffusive. However, if the perturbation in quark densities has a 
net non-zero electric charge density, then electromagnetic interactions can only be neglected 
if the scale of the fluctuation is small compared to the the electromagnetic transport mean 
free path of order (e 4 Tlne _1 ) _1 . On longer scales, the net charge density will relax at a 
rate affected by the electrical conductivity, while electrically neutral flavor asymmetries will 
relax diffusively. This will be discussed further momentarily. 

The leading-log diffusion constant characterizing the relaxation of fluctuations in charged 
lepton densities which are electrically neutral [e.g., an excess of electron minus positron 
density, balanced by an equal excess in p + — pT number density] is given by the appropriate 
specialization of Eq. ( 1.23 ) to QED, namely 



D L ~6 5 C(3) 2 7r- 



37T 5 

+ 4iV species + — 



e 4 In e 



;i.26) 



Once again, iV S pe C ies is the sum over all relevant fermion fields weighted by the square of 
their electric charge. In Eq. ( |L2£ ), and in the following Eqs. ( |1.29| )- (|1.31| ), the first term 



in the square bracket arises from t channel scattering from a gauge boson, the middle term 
represents t channel scattering from something else, and the third term arises from Compton 
scattering and annihilation to gauge bosons. 

The relaxation of an arbitrary set of slowly-varying fluctuations n a in net quark and 
(charged) lepton densities (with a labeling both quark and lepton species) is described by 
the coupled set of diffusion/relaxation equations 

^ = £> V 2 n --£e 6 n 6 , (1.27) 
ot e a Y 

where a a is the contribution to the conductivity due to charge carriers of species a (so 
that e a j a = cr a E, where j a is the species a particle number flux, and the total conductivity 
° — Ha a a)- These equations encode the fact that, in addition to various diffusive processes, 
the charge density p = J2a e a n a satisfies the non-diffusive relaxation equation p = —up, up 
to second order gradient corrections, showing that the conductivity is the relaxation rate 
for large-scale charge density fluctuations.^ In fact (as noted by Einstein), the conductivity 
is directly related to the underlying diffusion constants of individual species through the 
simple relation^ 



10 This, of course, immediately follows from combining the continuity equation for electric charge p+V-j = 0, 
the defining relation for conductivity j = crE, and Gauss' law V • E = p. 



11 This identity may be seen directly from the Kubo relations ( 1.8c ) and 
physical derivation (specializing for convenience to electromagnetism with a single species) is easily given. 



1.8d). Alternatively a simple 



Start with the diffusion equation j = —eD"Vn 
field the effective chemical potential is fi = (Mi 



-eD(dn/dfi)'V fi. 
eE • x. Hence j 



Then realize that in a constant electric 
2 D(dn/dfj,)E : or a = e 2 D{dn/dp). 
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(1.28) 



[or o a = e 2 a D a (dn a /dfi a )]. This relation assumes that the the charge susceptibility matrix 
( |1.9| ) is diagonal, as it is when all charge densities vanish. For (effectively) massless fermions, 
dn/d/j = \ T 2 per Dirac fermion. Since D L ^> D F , leptons completely dominate over quarks 
in the above species sum. Inserting the result ( |1.26| ) for the lepton diffusion constant into 
the Einstein relation ( |1.28|) reproduces our previous expression (|1.11|) for conductivity, as it 
must. 

For temperatures large compared to My/ (i.e., in the high temperature electroweak 
phase), one may find corresponding results for lepton number diffusion by specializing the 
general result (|1.23| ). If Yukawa interactions are neglected, then the relaxation of left and 
right-handed lepton number excesses are independent. The diffusion constant for left-handed 
net lepton number in high temperature electroweak theory is controlled by the SU (2) l gauge 
interactions and (approximately) equals 



6 5 C(3) 2 tt- 



6 + - (Nf + 2n s ) + — 
4 V ; 128 



1.29) 



with relative corrections of order l/lng" 1 , tern 2 9\y = {g'/dw) 2 , and a B . Here, N { = 12 
denotes the number of SU (2) chiral doublets, and n s is the number of scalar doublets. The 
corresponding diffusion constant for right-handed lepton number depends on hypercharge 
interactions, 



^ Lin — 



6 5 C(3) 2 tt- 



+ (20 + 2n s ) + 



3tt 2 



-i 



T- 



i-i 



1.30) 



with relative corrections of order l/lng' -1 and a s . Inclusion of Yukawa interactions will 
cause the diffusion of right and left-handed lepton number excesses to become coupled. 
This, however, only becomes relevant on scales larger than the mean free path for scattering 
processes involving Higgs emission, absorption or exchange. In the minimal standard model, 
this scale is of order [(mi/M w ) 2 (m t /M w ) 2 a 2 v T}~ 1 , where m e is the (zero temperature) mass 
of the lepton species of interest, and m t is the top quark mass. In other words, this scale is 
larger than the SU(2)l transport mean free path by a factor of roughly (M^/m^) 2 . 

The baryon diffusion constant for T ^> My/ is still given by the previous result ( |1.22| ) 
[with g = g s ], up to relative corrections of order 1/ In g" 1 . This may be rewritten as 



D 



B 



6 5 C(3) 



2 ^-3 
7T 



16 



16 Ar 2tt 2 

— N g H 

3 9 3 



-i 



T 



-i 



;i.3i) 



where N g = 3 is the number of generations. 

These (approximations to leading-log) diffusion constants for baryon and left or right- 
handed lepton number density characterize the relaxation in the high temperature elec- 
troweak phase of arbitrary fluctuations in any of these densities which are hypercharge 
neutral. For fluctuations having non-zero hypercharge density, one must also include the 
effect of the induced hypercharge electric field, leading to the same coupled relaxation equa- 



tions as in (|1.27|) , but with o now the hypercharge conductivity (and the species indices a, b 



12 



now labeling right-handed leptons, left-handed charged leptons, and quark flavors). For the 
hypercharge conductivity, the dominant contribution to the Einstein relation ( |1.28|) comes 
from right-handed leptons (since g 1 is the weakest coupling). Once again, one may easily 
check that inserting the result (|1.30| ) for the right-handed lepton diffusion constant into the 
Einstein relation reproduces our previous expression ( |1.13| ) for the hypercharge conductivity. 

We have not computed diffusion constants for conserved numbers carried by scalars. 

Several determinations of leading-log diffusion constants have previously been reported 
19| , |21| -p3| . Of these, only Moore and Prokopec |22 included all relevant diagrams, and 
each of these previous calculations made errors in evaluating at least one diagram. We will 
discuss the differences between our treatment and these previous results in more detail in 
Sec. 0. 



4. U(N f ) x U(N f ) flavor diffusion 

In QCD (or any QCD-like theory) with iVf species of fermions, one may consider the 
entire set of SU (N{ )y x SU (Nf )ax[/(1) b x £^(1)a currents. The diagonal components of the 
SU(N{)y currents are exactly conserved (neglecting weak interactions), while in our high 
temperature regime, the off-diagonal SU (iVf ) v currents are approximately conserved if one 
neglects order (Am/T) 2 effects, where Am is some fermion mass difference. Similarly, con- 
servation of the SU(N{) A currents is spoiled only by 0(m 2 /T 2 ) corrections. Consequently, 
the constitutive relations for the different currents must decouple, 

(1.32) 
(1.33) 
(1.34) 
(1.35) 

and the full diffusion constant matrix (in this basis) is diagonal when power corrections 
vanishing like T~ 2 , as well as weak and electromagnetic interactions, are neglected. 

The U(1)b baryon number current is almost exactly conserved,Q so fluctuations in baryon 
number density will behave diffusively on length and time scales large compared to the 
appropriate mean free scattering time. Fluctuations in flavor asymmetries — that is, the 
diagonal components of the SU(N{ )v current densities — will behave diffusively on time 
scales large compared to QCD mean free scattering times but small compared to the mean 
free time for flavor- changing weak interactions, which is of order (m\y /T) 4 / (a^T) , or (if the 
fluctuation is electrically charged) the electromagnetic transport mean free time of order 
(e^lne- 1 )- 1 . 

Fluctuations in the SU(N{) A and off-diagonal SU(Nf)y charge densities will behave 
diffusively on time scales large compared to mean free transport scattering times but small 



U(1) B : 


(JB> = 


-D B V(n B ), 


SU(Nt) w : 


Civ> = 


-D v V(n v ) , 


SU(Nf) A : 


0a) = 


-£>aV«>, 


U(1) A : 


(JA> = 


-D' A V(n A ), 



12 Baryon number is exactly conserved in QCD, but its conservation is violated by electroweak effects, at 
rates that are exponentially small at T <, Myy [p8 29 and 0(a^T In a" 1 ) at very high temperatures J30|-|32||. 
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compared to the time scale of order T/m 2 or T / {8m) 2 where the respective symmetry break- 
ing interactions become relevant.^ Because of the axial anomaly, fluctuations in the £7(1)a 
axial charge density may relax locally on a time scale of order (a^Tlna^ 1 ) -1 even in the 
massless theory. The physics behind this is completely analogous to treatment of baryon 
violating transitions in high temperature electroweak theory [|30|-|34|| .p^ As for the other 
axial currents, fluctuations in U(l) axial charge density will relax diffusively on time scales 
large compared to the QCD transport mean free times, but small compared to both the per- 
turbative T/m 2 and non-perturbative (af Tlna^ 1 ) -1 scales where {7(1)a violation becomes 
apparent. Axial charge fluctuations within this domain may be characterized by the basic 



diffusion equation (|1.35 ), with a diffusion constant D' A which is perturbatively computable. 

As will be discussed in detail in section |V|, to leading order in a s the various flavor 
diffusion constants only depend on two-to-two particle scattering rates in the high temper- 
ature plasma. Consequently, the diffusion constants for currents corresponding the various 
(approximate) flavor symmetry groups are all identical, 

D B = D V = D A = D' A , (1.36) 

up to relative corrections suppressed by one or more powers of a s . 



II. KINETIC THEORY AND TRANSPORT COEFFICIENTS 

To calculate any of the transport coefficients under consideration, correct to leading 
order in the interaction strength g but valid to all orders in 1/lng" 1 , it is sufficient to use 
a kinetic theory description for the relevant degrees of freedom. One introduces a particle 
distribution function /(p, x, t) characterizing the phase space density of particles (which one 
should think of as coarse-grained on a scale large compared to 1/T, but small compared to 
mean free paths). The distribution function /(p, x, t) is really a mult i- component vector 
with one component for each relevant particle species (quark, gluon, etc.), but this will 
not be indicated explicitly until it becomes necessary. The distribution function satisfies a 
Boltzmann equation of the usual form 



d d d 

777 Vp " — -|- I ox t ' 77 

at ax op 



f(p,x,t) = -C[f}. (2.1) 



The external force F ext term will only be relevant in discussing the electrical conductiv- 
ity. Since typical excitations in the plasma [those with p = 0(T)] are highly relativistic, 
corrections to their dispersion relations are suppressed by 0(g 2 ), and may be neglected. 



13 So in the presence of non-zero fermion masses, or mass differences, at sufficiently high temperature 
fluctuations in the SU(N{)y and SU(N{)a charge densities are "pseudo" -diffusive modes, analogous to 
pseudo-Goldstone bosons. Large scale fluctuations in these "almost-conserved" charge densities will satisfy 
a diffusion/relaxation equation of the form h = D\7 2 n — jn, where the local relaxation rate 7 will be 
0[(5m) 2 /T] or 0(m 2 /T), respectively. 

14 A useful discussion of this material may be found in Ref. Bq| [which, however, predates the realization 
that the transition rate per unit volume scales as 0(a 5 T 4 lna ), not as (aT) 4 ]. 
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Consequently, one may treat all excitations as moving at the speed of light, which means 
that the spatial velocity is a unit vector, v p = p = p/|p|-Q 

For calculations to leading order in g, and for the transport coefficients under consid- 
eration, it will be sufficient to include in the collision term C[f] only two-body scattering 
processes, so that 

C[f](p) = U \M(p,k,p',k')\ 2 (2 7 r)U 4 {p + k-p'-k') 

Jk,p',k' 



X 



{/(p) /(k) [l±/(p')] [1±/(10] - /(PO Z(k') [1±/(P)] [l±/(k)]} . (2.2) 



Here, p, k, etc., denote on-shell four-vectors (so that p° = |p|, etc.), M(p, k,p', k') is the two 
body scattering amplitude with non-relativistic normalization, related to the usual relativis- 
ts amplitude M. by 

|AW.*')r= , '^f.Xr (2 - 3) 

(2p )(2A;o)(2p , o)(2^) 

and / is shorthand for / c/ 3 p/(2-7r) 3 . The collision term is local in spacetime, and all dis- 
tribution functions are to be evaluated at the same spacetime position (whose coordinates 
have been suppressed). With multiple species of excitations there will, of course, be species- 
specific scattering amplitudes and multiple sums over species. As always, in the 1±/ final 
state statistical factors, the upper sign applies to bosons and the lower to fermions. Appro- 
priate approximations for the scattering amplitudes will be discussed in the next section. 
The stress-energy tensor, in this kinetic theory, equals 



T^(x)= / v£p»f(p,x), (2.4) 



p 



where v£ = p^/p° is a convenient generalization of the three- vector velocity for an excitation 
with spatial momentum p. (v£ is not the four-velocity and transforms non-covariantly, in 
just the manner required so that (d 3 p) v£ does transform covariantly.) Other conserved 
currents are given by similar integrals over the phase-space distribution function, but with 
the implied species sum weighted by appropriate charges q of each species,R 



f(x)= / v$qf(p,x). (2.5) 



p 



15 This assumes, of course, that the particular physical quantities under consideration are dominantly sen- 
sitive to the behavior of typical "hard" excitations. We will see that this is true for the transport coefficients 
under discussion. However, certain other observables, such as the bulk viscosity, may be sufficiently sensi- 
tive to the dynamics of "soft" excitations with momenta p <C T that their calculation requires an improved 
treatment which adequately describes both hard and soft degrees of freedom. 



16 Eq. (2.5) is adequate for currents which are diagonal in the basis of species. More generally, the distri- 
bution function /(p, x) should be viewed as a quantum density matrix for all internal degrees of freedom 
of an excitation. For example, in a theory with particles transforming in some representation R of the 
global symmetry group, the distribution function transforms under the R x R representation. Off-diagonal 
components of the distribution function are relev ant if one is interested in, for example, the off-diagonal 
parts of the SU{N { ) x SU(N f ) currents. Eq. (^|) would then be generalized to j£(x) = f v£ tr [i Q /(p, x)}, 
where t a is the appropriate charge representation matrix. 
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The factor p v which appears with /(p, x) in Eq. ( |2.4j) reflects the fact that in this case it is the 
energy or momentum of an excitation which is the conserved charge. Given the Boltzmann 
equation ( |2.1| ), and scattering amplitudes in ( |2.2| ) which respect the microscopic conservation 
laws, one may easily check that the currents (|2.4j) and Q2.5|) are, in fact, conserved. 

To extract transport coefficients, it is sufficient to linearize the Boltzmann equation (p.lfj 
and examine the response of infinitesimal fluctuations in various symmetry channels. This 
will be described explicitly below. 

But first we digress to discuss the validity of this kinetic theory approach. The Boltzmann 
equation fcj.lp may be regarded as an effective theory, produced by integrating out (off-shell) 
quantum fluctuations, which is appropriate for describing the dynamics of excitations on 
scales large compared to 1/T, which is the size of the typical de Broglie wavelength of an 
excitation. The use of kinetic theory for calculating transport coefficients may be justified 
in at least three different ways: 

1. One may begin with the full hierarchy of Schwinger- Dyson equations for (gauge- 
invariant) correlation functions in a weakly non-equilibrium state in the underlying 
quantum field theory. For weak coupling, one may systematically justify, and then in- 
sert, a quasi-particle approximation for the spectral densities of the basic propagators, 
perform a suitable gradient expansion and Wigner transform, and formally derive the 
above kinetic theory. (See Refs. [BB-RBI and references therein.) 



2. One may consider the diagrammatic expansion for the equilibrium correlator appearing 



in the Kubo relation (jl.8|) for some particular transport coefficient. After carefully 
analyzing the contribution of arbitrary diagrams in the kinematic limit of interest 
(k = and uj — » 0), one may identify, and resum, the infinite series of diagrams which 
contribute to the leading-order result. One obtains a linear integral equation, which 
will coincide exactly with the result obtained from linearizing the appropriate kinetic 
theory. This program has been carried out explicitly for scalar theories but not 
yet for gauge theories. 

3. One may directly argue (by examining equilibrium finite temperature correlators) that, 
for sufficiently weak coupling, the underlying high temperature quantum field theory 
has well-defined quasi-particles, that these quasi-particles are weakly interacting with 
a mean free time large compared to the actual duration of an individual collision, 
and consequently that scattering amplitudes of these quasi-particles are well-defined 
to within a precision of order of the ratio of these scales. In other words, one justifies 
the existence of quasi-particles by looking at the spectral densities of the propagators 
of the basic fields, reads off their scattering amplitudes from looking at higher point 
correlators, and writes down the kinetic theory which correctly describes the resulting 
quasi-particle interactions. For a more detailed discussion of this approach see, for 
example, Ref. pBJ. 



Given the complexities of real-time, finite-temperature diagrammatic analysis in gauge the- 
ories (especially non-Abelian theories), we find the last approach to be the most physically 
transparent and compelling. But this is clearly a matter of taste. 

There is one important caveat in the claim that a kinetic theory of the form ( |2.1| ) can 
accurately describe excitations in a hot gauge theory. One may argue, as just sketched, 
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that such a Boltzmann equation with massless dispersion relations reproduces (to within 
errors suppressed by powers of g) the dynamics of typical excitations in the plasma, namely 
hard excitations whose momenta are of order T. For such excitations, thermal corrections 
to the massless lowest-order dispersion relations are a negligible 0(g 2 ) effect. This is not 
true for soft excitations with momenta of order gT, or less. In gauge theories, one cannot 
characterize sufficiently long wavelength dynamics in terms of (quasi)-particle excitations 
with purely local collisions. Instead, one may think of long wavelength degrees of freedom 
as classical gauge field fluctuations, and construct Boltzmann- Vlasov type effective theories 
which describe hard excitations propagating in a slowly varying classical background field. 



The well known hard-thermal-loop (HTL) effective theory is of precisely this form [p9| - p2| , |36 

A simple kinetic theory of the form (|2.1|) , without the complications of background 
gauge field fluctuations, can only be adequate for computing physical quantities which are 
not dominantly sensitive to soft excitations. This is true of most observables, including 
thermodynamic quantities such as energy density or entropy, just because phase space grows 
as p 2 <i|p| in (3+1) dimensional theories. This is equally true for the transport coefficients 
under consideration. It will be easiest to demonstrate this a-posteriori. However, this 
insensitivity (at leading order) to soft excitations may not hold for the bulk viscosity, which 
is why its calculation requires a more refined analysis. (This is true even in a pure scalar 
theory @,|.) 

Returning to the analysis of the Boltzmann equation (|2.1| ), equilibrium solutions are 
given by 



/ e » = {exp \p{-u v tf - fi a q a a )\ T l}" 1 , (2.6) 

where (3 is the inverse temperature, u is the fluid four-velocity, and are chemical 

potentials corresponding to a mutually commuting set of conserved charges. We have now 
included an explicit species index a, and is the value of the a'th conserved charge carried 
by species a. (Sums over repeated charge indices should be tacitly understood.) 

Using only the fact that the scattering amplitudes respect the microscopic conservation 
laws, one may easily show that the collision term exactly vanishes for any such equilibrium 
distribution, C[/ eq ] = 0. 

The distribution function corresponding to some non-equilibrium state which describes 
a small departure from equilibrium may be written as the sum of a local equilibrium distri- 
bution plus a departure from local equilibrium. This is conveniently written in the form 

HP, x) = / a (p, x) + / a (p, x) [1 ± / a (p, x)] #(p, x) . (2.7) 

Here Jq(p,x) has the form of an equilibrium distribution function, but with temperature, 
flow velocity, and chemical potentials which may vary in spacetime, 

Writing the departure from local equilibrium as /o(l±/b) fi, instead of just 5f, simplifies the 
form of the resulting linearized collision operator ( |2.10|) . When inserted into the collision 
term of the Boltzmann equation, the local equilibrium part of the distribution gives no 
contribution, C[fo] = 0, because the collision term is local in spacetime and so cannot 
distinguish local equilibrium from genuine equilibrium. Hence, the collision term, to first 
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FIG. 1. Momentum and species label conventions for 2 — > 2 scattering. Time runs from left to right. 

order in the departure from equilibrium, becomes a linear operator acting on the departure 
from local equilibrium, 



(2.9) 



C\f\=Cf 1 + 0{ft), 
where the action of the linearized collision operator C is given by 

ffhc . „ 

W(p) = |E M%(p,k,p',k') [2^5\p + k - jf - k') 

bed Jk 'P'' U ' 

x fS (P) /o 6 (k) [1±/ C (P')] [l±/o(k')] [/f(p) + A b (k) - /r(p') - /f(k')] . (2.10) 

All distribution functions are evaluated at the same point in spacetime, whose coordinates 
have been suppressed. M^(p, k,p', k') is the scattering amplitude for species a and b, with 
momenta p and k, respectively, to scatter into species c and d with momenta p' and k'. 
For reference, this choice of momentum and species labels is summarized in Fig. [l|. The "c" 
in the "ffhc" above the sum indicates that, in the application to gauge theories, the sum 
is over all colors of the particles represented by b, c, and d, as well as flavors, anti-flavors 
(where distinct), and helicities. 

On the left-hand side of the Boltzmann equation, the gradients acting on f give a result 
whose size is set by either the magnitude of spacetime gradients in temperature, velocity, or 
chemical potentials, or by the magnitude of the imposed external force. So, to first order in 
the departure from equilibrium, the Boltzmann equation becomes an inhomogeneous linear 
integral equation for /i(p), 



<9 „ <9 d 
dt +P 'd^ + cxt 'd^ 



/ a (p,x,t) = -(CAr( P ,x,t) 



(2.11) 



In other words, fi is first order in gradients (or the external force). The neglected terms 
on the left-hand side, where derivatives act on the deviation from local equilibrium, are of 
second order in gradients (or external force), and do not contribute to the linearized analysis. 

Each transport coefficient under consideration will depend on the departure from equi- 
librium fx resulting from a particular form of the driving terms on the left-hand side of the 
linearized Boltzmann equation fl2.11|) . For conductivity, one is interested in the response 
to a homogeneous electric field, where Fg xt = q a ~E. For diffusion or shear viscosity, one is 
interested in the response to a spatial variation in a chemical potential or the fluid flow 
velocity. Using the fact that 



d/ (p, x) = /o(p, x)[l ± / (p, a;)] d[0u u p v + /3fj, a q a a ] 



(2.12) 
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one may easily see that in all three cases, the left hand side of the linearized Boltzmann 
equation ( 2.11 ) has the form,^] 



LHS = PfS(p,x)[l ± / a (p,x)] ,," /,... / (p).V,.., ; (.r) , (2.13) 
where the spatial tensor Xi...j(x) denotes the "driving field," namely 

{ — Ei , (conductivity) 
Vi/i a , (diffusion) (2.14) 

^7= (ViUj + VjUi — §<5ijV ■ u) , (shear viscosity) 

and Ii...j(p) is the unique £ = 1 or £ = 2 rotationally covariant tensor depending only on the 
direction of p, that is 

pi, (conductivity/diffusion) (2 15) 

S (PiPj ~ ■ (shear viscosity) 



In Eq. (|2.13|) , q a denotes the relevant charge of species a which, in the case of shear viscosity, 



means the magnitude of its momentum |p|. The factor of w| included in the definition 



( |2.15| ) of Iij [and the corresponding 4g factor for in Eq. (|2.14|) 1 are inserted so that the 



normalization 1^1^ = 1 holds for both the £ = 1 and £ = 2 cases. 

We will henceforth always work in the local fluid rest-frame. At any point x, the local 
equilibrium distribution function /o(p, x) then depends only on the energy p° = |p|, and 
thus the only angular dependence on p in ( [2.13] ) comes from the £ = 1 or £ = 2 irreducible 
tensor Jj...j(p). Since the linearized collision operator C is local in spacetime and rotationally 
invariant (in the local fluid rest-frame at x), the departure from equilibrium must have the 
same angular dependence as the driving term. Consequently, given a left-hand side of the 
form (|2.13|) , the function /i(p,£) which will solve the linearized Boltzmann equation (\2.11 



must have the corresponding form 

A a (p,x) =p 2 X l ... J (x) /^-(p) X a (|p|), (2.16) 

where, for each species a, X a (|p|) is some rotationally invariant function depending only 
on the energy of the excitation. The factor of (3 2 is inserted for later convenience, and 
causes x a to have the same dimensions as q a (dimensionless for conductivity or diffusion, 
and dimension one for viscosity). For notational convenience we will also define 

Xt.,(p) = /,..,(p)x a (|p|), (2.17) 

and 

^■(p) = -Tq a fg(p,x)[l ± / a (p,x)] J,..,(p) , (2.18) 



1 For diffusion or viscosity, the time derivative term on the left-hand side may be dropped because its 
contribution is actually second order in s pati al gradients. This follows from current (or stress-energy) 
conservation and the constitutive relations ([0|) or ([I]?]), which together imply that time derivatives of the 
conserved densities are related to second spatial derivatives of the conserved densities themselves. 
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so that the linearized Boltzmann equation for any particular channel under consideration 
can be written in the concise form 



st.Jp) = (Cxi-iT (p) ■ 



(2.19) 



A straightforward approach for numerically solving these coupled integral equations 
would be to reduce them to a set of scalar equations [by contracting both sides with Jj...j(p)], 
discretize the allowed values of |p|, compute the matrix elements C ab (|p|, |q|) of the kernel of 
the (projected) collision operator by numerical quadrature, and thereby convert ( |2.19|) into 
a finite dimensional linear matrix equation. This is a bad strategy, however, particularly for 
gauge theories. The problem is that the kernel C ab (|p|, |q|) has (integrable) singularities and 
is not smooth as |p| crosses |q|. Consequently it is quite difficult to avoid large discretization 
errors and obtain good convergence to the correct answer. 

A better strategy, which is applicable to the full leading-order analysis and not just 
the present leading- log treatment, is to convert the linear integral equations ( |2.19| ) into an 
equivalent variational problem. This permits one to obtain quite accurate results using very 
modest basis sets. This conversion is trivial once one notes that the linear operator C is 
Hermitian with respect to the natural inner product 



ffhc „ 

(f, 9 ) = /3 3 Y: J 



Consequently, if one defines the functional 

Q[x] = (xi-j,Si...j) - ~ (xi-j,CXi-j): 



(2.20) 



(2.21) 



then the maximum value of Q[x] occurs when x a (|p|) satisfies the desired linear equation 
( |2.19j ). For later use, note that the maximal value of the functional Q may be written in 
either of the forms 



X=Xn 



2 (Xi-j) 



Qma,x 2 \X.i--j:^X 

In more explicit form, the two terms in Q are 

ffhc 

(x*..i, #...;) = -/3 2 £ / / (p)[l±/o(p)] g°X a (|pi), 

VXi—j : ^Xi-'-j 



X=Xn 



[2.22) 



(2.23) 



and 



53 ffhc 



8 i^P-k-P'-k' 



M%(p, k, p', k') (2vr) 4 5\p + k-p'-k') 



x fS(p) / b (k) [l±/ c (p')] [l±/ d (k')] 



X 



X-.,(P) + Xt,(k)-X-.,(P')-Xt,(k') 



(2.24) 



The sum is over all scattering processes in the plasma taking species a and b into species 
c and d. We have used crossing symmetry of the scattering amplitudes to write the above 
expression for matrix elements of C in a form which makes it apparent that C is a positive 
semi-definite operator. The overall factor of 1/8 compensates for the eight times a given 
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process appears in the multiple sum over species when all the particles are distinct (due to 
relabeling a b, c d, and/or ab <-> cd), and supplies the appropriate symmetry factor 
in cases where some or all of the particles are identical. For our later discussion it will 
be important to note that C is non-diagonal in the basis of species when there are 2 <-> 2 
processes involving more than one species type. 

After solving the linearized Boltzmann equation in the particular channel of interest, 
by maximizing Q[x], the associated transport coefficient may be determined by inserting 
the resulting distribution function [given by Eqs. ( |2.7| ) and (|2.16|) 1 into the stress-energy 



tensor fl2.4p or appropriate conserved current ( |2.5| ), and comparing with the corresponding 
constitutive relation [Eq. (|1.4|), (|1.6|), or (|1.7|) 1- In each case, the integral over the distribution 



function which defines the flux (i.e., the stress tensor Tjj or a spatial current jf) reduces to 
the inner product (xi-j,Si...j). Consequently, the actual value of each transport coefficient 
turns out to be trivially related to the maximum value (|2.22| ) of the functional Q in the 
corresponding channel. Explicitly, 

(2.25) 
(2.26) 

= 1, q=q a yOfiaJ 

The thermodynamic derivative appearing in fl2.26| ) is the charge susceptibility, 

dn r fihc 




d/j,, 



- = E(^) 2 //?/o(p)[1±/o(p)] 

a a J P 



ffhc 



^ 2 £A a (g*) 2 , (2.28) 



where, once again, A a is 1 for fermions and 2 for bosons, and we have specialized to the 
ultra-relativistic limit. 



III. COLLISION INTEGRALS 

The full set of scattering processes which contribute at leading order in a theory with 
gauge and fermionic degrees of freedom are shown in Fig. [|. Some of these processes yield 
matrix elements which become singular as the momentum transfer (i.e., Mandelstam t or 
u) goes to zero. For instance, in a vector-like theory, the matrix element for gauge boson 
exchange between fermions [diagram (C)} behaves, for non-identical fermions, as 

which diverges as the inverse fourth power of the momentum transfer as p approaches p'. 
Phase space only partially compensates, resulting in a cross section which is quadratically 
divergent at small p — p'. Similarly, the annihilation diagram (D) has Ai 2 oc (u/t) + (t/u), 
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FIG. 2. Leading-order diagrams for all 2 <-> 2 particle scattering processes in a gauge theory with 
fcrmions. Solid lines denote fermions and wiggly lines are gauge bosons. Time may be regarded as running 
horizontally, either way, and so a diagram such as (D) represents both // — > gg and gg — > //. The diagrams 
of the first row [(A)— (E)] contribute to the leading log transport coefficients, while the diagrams of the second 
row [(F)-(J)], and all interference terms, do not. 

which leads to a logarithmically IR divergent scattering cross section. However, we are not 
directly interested in the total scattering cross section; we need to know the size of the 
contribution to (Xi-jiCXi-j) i n the channels relevant to transport coefficients. As we shall 
review, these transport collision integrals can be less singular than the total scattering rate. 



A. Kinematics 



It is convenient to arrange the phase space integrations so that the transfer momentum is 
explicitly an integration variable. This will make it easy to isolate the contribution from the 
potentially dangerous small momentum exchange region. We choose to label the outgoing 
particles so that any infrared singularity in (a given term of) the square of the amplitude 
\A4\ 2 occurs only when (p'—p) 2 — > 0.0 In the collision integral ( (2.24|) it is convenient to 
use the spatial 5 function to perform the k' integration, and to shift the p' integration into 
an integration over p'—p = q. We may write the angular integrals in spherical coordinates 
with q as the z axis and choose the x axis so p lies in the x-z plane. This yields 



\Xi-j-i CXi—j) 



Ij3 ffhc 

(4tt) 6 



ffhc poo rl p2TT 

/ q 2 dq p 2 dp k 2 dk \ d cos 6 vq d cos 6k q \ 
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5( P +k-p'-k') fs(p) m [i ± / v)] [i ± f d (k" 



X 



l-.J 



:P) + xt,(k)-X-.,(p / )-xt ? (k / ) 



(3.2) 



18 There is one case where this is impossible, namely, scattering between identical fermions, where the 
interference term between outgoing leg assignments in diagram (C) makes a contribution to the matrix 
element M 2 oc s 2 /ut, which is divergent for both t — > and u — * 0. As will be discussed shortly, this 
interference does not contribute at leading-log level. Regardless, one could also put this case in the desired 
form by using s = —u — t and rewriting the matrix element (squared) as (s/t) + (s/u), so that each piece is 
now singular in only one momentum region. Diagram (A) apparently has the same problem; but when one 
sums all gg — > gg processes (only the sum is gauge invariant) one finds M 2 oc (3 — su/t 2 — st/u 2 — tu/s 2 ), 
so there is no problem. 
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where here and henceforth, p, k, and q denote the magnitudes of the corresponding three- 
momenta (not the associated 4-momenta), p' = |q + p| and k' = |k — q| are the magnitudes 
of the outgoing momenta, <fi is the azimuthal angle of k (and k') [i.e., the angle between 
the p-q plane and the k-q plane], and 9 pq is the plasma frame angle between p and q, 
cos 9 pq = p • q, etc. 

Following Baym et al. [ 17fl , it is convenient to introduce a dummy integration variable 
u, constrained by a 5 function to equal the energy transfer p' — p, so that 



/oo 
duo 5(lu + p — p) 5(lo — k + k!) . 
-oo 



Evaluating p' = |p + q| in terms of p, q, and cos 9 pq , and defining t = 
usual Mandelstam variable), one finds 

8(00 +p-p') = —d(cos9 pq - — - — ) <d(uj + p) , 
pq V q 2pq) 

kl / uo t \ 

5(uo-k + k') = — 5[cos9 kq - - + —J 0(k - uo) , 



(3.3) 

uo 2 — q 2 (which is the 

(3.4) 
(3.5) 



where 0(z) is the unit step function. The cos 9 integrals may now be trivially performed 
and yield one provided p > |(g — uj), k > |(g + u), and \u\ < q; otherwise the argument of 
a S function has no zero for any | cos#| < 1. The remaining integrals are 
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with p' = p + u, and k' = k — uo. For evaluating the final factor of 

Ii...j(i>) 4..j(k) = Pe(cos9 pk ) , 



note that 



(3.6) 



(3.7) 



where Pi is the £'th Legendre polynomial. We will therefore need expressions for the angles 
between the momenta of all species, as well as the remaining Mandelstam variables s and 
u, which appear in M 2 . They are 
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(3.8) 
(3.9) 

(3.10) 

(3.11) 
(3.12) 
(3.13) 
(3.14) 
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and 



s = 2pk(l-cos6 pk ) = — {(p+p')(k+k') + q 2 - cos -/(4V + t) (4fcfc' + t) j , (3.15) 
u=-t-s. (3.16) 



B. Leading-log simplifications 

To compute leading-log transport coefficients, it will be sufficient to extract the small q 
contribution to the collision integral ( |3.6| ). For small q but generic p and k (so that q <T, 
q ^ p, and q -C k), one has 

— s u 2pk .„ . „„. 

— ~ - ~ -cos0). (3.17) 

t t q 2 y 1 y 1 

The integral is restricted to the range —q<uj<q and will be dominated by cu/g of order 
one. Hence, the small q integration region in a diagram where M 2 oc (s 2 or u 2 )/t 2 naively 
behaves like / dq/q 3 , a quadratic infrared divergence; while for a diagram with Ai 2 oc st/t 2 
(or wt/t 2 or s 2 /tu) the small g integration region naively behaves like J <ig/g, a logarithmic 
divergence. These estimates are inadequate, however, if species a and c are identical, and 
species b and d are also identical — that is, when both incident particles undergo small angle 
scattering without changing their species types. In this case, since p' — p = q is small, one 
has 

Xl,(p) - Xtjip') = -q ■ Vxt.;(p) + 0(q 2 ) , (3.18) 

and similarly for Xi-jO*) ~ Xi-jO^') ■ Therefore, the [x a +X b ~ X°~ X d } 2 factor in the collision 
integral (|3.6j ) will contribute a factor of q 2 to the integrand, softening the small q divergence. 

This q — > cancellation is operative for diagrams (A), (B), and (C), and converts the 
naive estimate of a quadratic infrared divergence into a merely logarithmic divergence. It 
also converts interference terms involving these diagrams, which were naively log divergent, 
into finite contributions. The cancellation does not occur for diagrams (D) and (E) which 
involve a change of species. For these diagrams, the naive estimate of an infrared log 
divergent result is correct. Interference terms involving these diagrams, and the remaining 
diagrams (F)-(J), are infrared finite from the outset. Hence, (the squares of) all diagrams 
in the first row in Fig. ^| lead to logarithmic IR divergences in Eq. ( p.6| ), while diagrams in 
the second row, and all interference terms, do not. 

All of these logarithmic divergences become convergent when one includes the self- 
energies appearing on exchange lines in the diagrams. The self-energies are all of order 
g 2 T 2 . For small 5 < T, the g 2 T 2 part of the self-energies are known as the hard thermal 
loop (HTL) self-energies [41 , 42| . An 0(g 2 T 2 ) self-energy correction on a propagator is im- 



portant when the exchange momentum squared becomes of order g 2 T 2 , which means when 
q ~ gT. In every case, the inclusion of the self-energy reduces the size of the matrix element 
and serves to cut off the log divergence in the infrared. For the case of the gauge boson 
propagator, relevant in diagrams (A)-(C), this is discussed in [0. The demonstration that 
the log is cut off for diagrams (D) and (E) has apparently not appeared in the previous 
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\i,a 

FIG. 3. A generic vertex from diagrams (A)-(C) of Fig. ^, to be analyzed in the soft exchange limit. In 
this figure, the solid line denotes any sort of particle (e.g., a gauge boson or fermion) that is being scattered, 
and the wavy line represents the exchanged gauge boson. 



literature, though the self-energy for the fermion line is well known |]43 |. Our analysis shows 
that the self-energy on the fermion line is sufficient to cut off this IR divergence as well ||. 
In the current paper we will not discuss this issue in detail; nor will we treat carefully the 
momentum region q > T, where the small q approximations which led to the conclusion that 
there is a log divergence break down, cutting off the log from above. Instead we will make 
a leading-log treatment, which means that we will extract the coefficient of the logarithmic 
divergence. This permits us to simultaneously take q to be small, q <C T, q <C p, q <C k, 
allowing certain kinematic simplifications, while simultaneously neglecting self-energy cor- 
rections in determining the matrix elements. This approximation has been customary in 
almost all work in the field; in fact we are not aware of any paper which correctly goes 
beyond this approximation when computing a transport coefficient in a relativistic gauge 
theory. In a companion paper, we will treat the problem to full leading order in g §. 

C. Diagrams (A), (B), and (C) 

Consider the gauge boson exchange diagrams (A), (£?), and (C). These represent oh <-> oh 
processes where, in the relevant small q regime, the incoming and outgoing lines with nearly 
the same momenta are the same species type. The near cancellation ( |3.18| ) thus implies that 
the x factors will contribute an explicit q 2 for small q. We may therefore use leading small 
q approximations everywhere else, together with the first nontrivial small q approximation 
for the Xt-j factors. As a result, the various forms which appear in the square of the matrix 
elements for the different diagrams (A)-(C) all become the same to leading order. In fact, 
in the soft exchange (q — > 0) limit, the vertices in diagrams (A)-(C) take on a universal 
form that depends only on the color charge of the particle that is scattering, be it a gluon or 
quark (or scalar). Such a vertex is depicted generically in Fig. [| and the associated vertex 
factor in the q — > limit is0 

2p»gt a 5 hh ,, (3.19) 

where t a is the color generator in the representation of the scattering particle, and h and h! 
are the ingoing and outgoing helicities of that particle. (The 2p in the factor 2p M = 2pv 



19 One may alternatively work directly with the full matrix elements and observe that their structure 
becomes identical in the q — > limit because (s 2 +u 2 )/(2t 2 ) ~ —su/t 2 ~ s 2 /t 2 ~ (Ap 2 k 2 /q i ) (1 — cos<ji>) 2 . 
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is a consequence of using relativistic normalization for matrix elements.) With this q — > 
Feynman rule, it is trivial to evaluate the matrix elements for the ab «-> ab scattering 
processes of Figs. (A)-(C): 



M 



leading- 



A a b 



(3.20) 



where the coefficient A a b depends on the gauge coupling and representations of the two 
species,^ 



Aab = 4:d A T Ra T Rb g 4 



(3.21) 



Strictly speaking, this gives the square of the amplitude summed over incoming and outgoing 
gauge group indices and outgoing spins or helicities, but not summed over incoming spins 
or helicities, which are considered part of the species label. The corresponding contribution 
to the collision integral becomes 
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x fS(p) [1 ± f a M\ foik) [1 ± f (k)} 



xl,(p) " X-.,(P) + UL 7 (k) - Xt,(k') • (3.22) 



x 



As in Eq. ( |L"23|) , the species sums run over all helicities and types of excitations, count- 
ing anti-particles separately from particles. The 2 next to A ao arises because the sum in 
Eq. ( [2.24|) separately counts both ab — ► ab and ab —> 

The leading small q approximation (|3.18 ) to the difference of \ factors, for either i = 1 
or £ = 2, has the explicit form 



x? (p) - xf(p') = -whip) x a {v)' + (wA - ft) 



x a (p) 



p 



+ 0(g 2 



Xy(p) - X?i(pO = -w/i/(p) X a (p)' + Ji (ZuiPiPj ~ qiPj - qjPi) 



x a (p) 



P 



(3.23a) 
+ 0(q 2 ) , (3.23b) 



where x a (p)' means dx a (p)/dp. Expressions for x|..,-(k) — xt-jO^) are ^ ne same except for 
replacing p by k, and changing the overall sign. For either case (and in fact, for any £), 

(3.24) 



xl,(p) - xL.Ap)} 2 = lx a ( P )f + q ^-^ ix a (p)] 2 + 0(g 3 ) • 



P 



20 If a=b, then the matrix element has a second s 2 /u 2 term arising from the interchange of the outgoing 
lines. Swapping the p' and k' labels effectively makes the resulting matrix element twice as large. 

Unless a=b, in which case the extra factor of two comes from the matrix element, as noted in the previous 
footnote. 
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When expanding the last factor of Eq. ( |3.22| ), there are two types of contributions to 
examine: those involving two \ a or two \ b factors, for which one may use the expression 
( p.24| ), and the cross-contributions with one x a an d one \ b ■ We will consider the cross- 
contributions first. As noted above, the explicit q or uj appearing in the difference (]3.23|) 
softens the small q behavior to a logarithmic divergence, so in all other factors one may work 
to leading order in q. In particular one may approximate 



OJ 

COS 6 pq ~ COS Okq — COS 9 p i q ~ COS Ok'q — — , 

uj 2 q 2 — ijj 2 

cos 6 p k — cos p 'k — cos 6 p k> — cos 6 p 'k' — — H ^ — cos 4> ■ 



q- 



(3.25) 
(3.26) 



Explicitly carrying out the du and dtp integrations, with the (1 — cos0) 2 factor from the 
matrix element included, one finds that all x a X b cross terms vanish in case of I = 2 (or 
higher), but not for £ = 1. In the t = 1 channels, however, we will only be interested in 
the diffusion of charge conjugation (C) odd quantum numbers (electric charge or various 
fermionic numbers like baryon number). In our high temperature regime, the equilibrium 
state may be regarded as C (or CP) invariant. Consequently, C (or CP) symmetry ensures 
that particles and anti-particles will have opposite departures from equilibrium, \ a = —x a 
(where a denotes the anti-particle of species a). Hence, the sign of the cross term will be 
different for scatterings from fermions versus anti-fermions, so that the two contributions will 
cancel in the sum over species. When the scattering involves a gauge boson on one or both 
lines, then C symmetry ensures that the departure from equilibrium for the gauge boson 
is zero, so again there is no cross term. For the same reason, gauge-boson — gauge-boson 
scattering [diagram (A)} plays no role for conductivity or diffusion. 

In either case, what remains are only the terms with two x a or two \ b factors. After 
using (|3.24j) (or the same relation with a «-> b and p <-> k), the u and <fi integrals are simple, 
and the k integral can also be performed using 



/ dkk 2 f^k)[l±f b (k)] = \ b T 
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(3.27) 



where Af, = 2 if species b is bosonic, and 1 if it is fermionic. The result is 

\ diagrams (A—C) 
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X b / » [1 ± r {p)} (p 2 X a {p)' + X a {p? ) + (a~b) 



(3.28) 



In the q integration, the upper cutoff is q ~ T, where the small q treatment breaks down. 
The lower cutoff occurs because we have not included the gauge boson's hard thermal loop 
self-energy in computing the matrix element. Inclusion of the self-energy makes the matrix 
element smaller than the vacuum amplitude and cuts off the q integration in the infrared.0 
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This is true for both longitudinal and transverse parts of the exchanged gauge boson propagator [fi"7f . 
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Hence, in a leading log treatment one may simply replace the entire q integral by log g , 
and thereby reduce these collision integral contributions to a single one- dimensional integral. 



D. Diagrams (D) and (E) 

We begin with the annihilation diagram (D). The matrix element squared for a fermion 
of species / to annihilate with an anti-fermion of opposite helicity and produce two gauge 
bosons, // — > gg, summed over initial and final gauge group indices and gauge boson spins, 
is 



M ff 



leading— log 
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u t 
t u 



(3.29) 



with 



A f =4d A T Rf C Rf g 4 



(3.30) 



Interchanging labels on the outgoing legs turns (t/u) into (u/t), and so one may keep just 
the u/t part of the matrix element and multiply the result by two [which effectively cancels 
part of the overall 1/8 symmetry factor which appears in Eq. ( |2.24j )1. 

Since the degree of divergence is at most logarithmic, we may immediately make all 
available small q approximations. Namely, we may take the limits of the p and k integrations 
to be zero, take f (p+u>) ~ /o(p) and similarly f Q (k—uj) ~ fo(k), approximate cos 6^/ ~ 
cos 9kk' — 1, an d use Eqs. ( |3.25| ) and ( |3.26| ) for the various angles. The matrix element, at 
leading order in small q, is just 



u 



t 



2kp 
q 2 



cos « 



(3.31) 



Making these approximations gives the following contribution to the collision integral ( |2.24| ) , 
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X f (k)-X 9 (k) 



X f (p)-X 9 (p)}[x f (k)-X 9 (k)}}- (3.32) 



Here x is the departure from equilibrium for the fermion, x is f° r its anti-particle, and x g 
is for the gauge boson. The distribution function fl is the equilibrium Fermi distribution 
while /q is the equilibrium Bose distribution. The sum runs over all fermion species and 
helicities, but not over anti-particles. The factor of 8 next to Af is the 2 from the two pieces 
of the matrix element, times the 4 ways Eq. Q2.24j ) counts this diagram (// — > gg, ff — > gg, 

gg -> ff, gg^Jf)- 

Using Eq. ( |3.26| ) for Q pk , one may easily check that 
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/q rZw 
duo J d(j) (1 — cos0) Pe(cos9 p k) = (3.33) 

for £ = 1 and £ = 2 (and in fact, all £ > 0). Therefore the cross term involving both x^ 
and vanishes. For the remaining terms, the u and <fi integrals are trivial. For the term 
involving [x^{k)—X g {k)] 2 , one may use the symmetry of the integration region to exchange 
p and k. After performing the k integration with 

| o °° dk k fi(k) [i + = ^ , (3.34) 

one finds 

V J / leading-log *-y 2 y 7T d J gT q JO 

x {[x f (p) - X 9 (P)] 2 + [X 7 (P) - X 9 (p)f) ■ (3-35) 

Once again, the limits on the q integration show where the approximations we have 
used break down. The small q approximation breaks down for q > T. Approximating the 
matrix element by the vacuum amplitude, without a thermal self-energy insertion on the 
internal fermion propagator, is invalid for q < gT. Below this scale the matrix element is 
smaller — parametrically smaller for q gT. Our approximations are an over-estimate for 
q < gT, and the integral is cut off. For a leading log treatment, one may simply replace 
the q integration with \og(T/gT) = \ogg~ 1 , which reduces the collision integral contribution 
( j3.35|) to a single integral over a quadratic form in the departures from equilibrium. 

Finally, the Compton scattering diagram (E) differs only slightly from the annihilation 
diagram (D). The matrix element for diagram (E) is —(s/t) rather than (u/t), but at 
leading order in small q these are equivalent. The sign of each x associated with excitations 
with momentum k is opposite, and the fermion departures from equilibrium are now either 
both x^ or both x i rather than one of each. Summing both particle and anti-particle 
cases, and exploiting the fact that the p-k cross-term cancels, one finds that the leading log 
contribution of diagram (E) is identical to that of diagram (D), 



( n \ diagram (E) / \ diagram (D) /oor\ 

[Xi-J, Cxi-j J leading „ log = (Xi-i, C ^...ij leadiiig _ Iog • (3.36) 



Compton scattering and annihilation/pair-creation processes do not give equal contributions 
beyond leading log, but that does not concern us in the present paper. 



IV. ELECTRICAL CONDUCTIVITY 



To find the electrical conductivity, one must determine the departure from equilibrium 
which is produced by an imposed electric field. For the DC conductivity, one may take the 
electric field E to be constant and neglect all spatial and temporal variation in the distribu- 
tion function. As discussed in section ||, only the external force F" xt = q a E contributes on 
the left-hand side of the linearized Boltzmann equation ( |2.11| ). Rotational invariance implies 
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that the correction /"(p, x) to local equilibrium may be taken to have the form (|2.16|) which, 
specialized to the current case, reads 



f?(p,x) = -p 2 (E-p) X a (p). (4.1) 

The linearized Boltzmann equation reduces to the form (p) = (Cx«) a (p) > with Sf(p) = 
-Tq a fS(p) t 1 ± /o (p)] Pi and Xi (p) = Pi X a (p)- The conductivity is given by 

• Tp -p ffhc „ 

o- = = ^ • £ g. / p p /o (p) [i - /o (p)] A a (p) 

ffhc „ 

= -|/5 2 E^/ /o(p) [l-/o(p)]x a (p) 



2 Qmax|^ = i, 0=0 „.. , (4.2) 



- 3 ^ max 1^=1 ,q=q EM , 

where Q max is the maximal value of the functional Q[x] = (Xi> <%) ~ | (x«> ^Xi), first defined 



in Eq. ( 2.21 ). [In the explicit forms shown in Eq. ( |4.2|) , we have specialized to only fermionic 
charge carriers, as is appropriate for both pure QED and the full standard model.] 

The driving term Si(p) in the linearized Boltzmann equation is charge conjugation (C) 
odd, and therefore the departure from equilibrium x(p) must likewise be C odd.0 Hence, 
\ e (p) = — x e+ (p), while x 7 (p) = 0- 111 our high temperature regime, the departure from 
equilibrium will also be identical for different helicities of the same species, and for different 
leptons of the same charge if there are multiple light lepton species, so that x M+ = X e+ ; e ^ c - 
If there are active quark species, then rapid QCD scattering processes (on the time scale of 
the relevant electromagnetic interactions) ensure that X q (p)/x e+ (p) = up to 0[(ckem/ 'ot s ) 2 ] 
corrections, which we will neglect. (Active quarks remain relevant, however, as targets from 
which the charged leptons can scatter.) Hence, for this application we may regard x e+ ip) as 
the only independent function which must be determined. 

Let iVieptons denote the number of active lepton species (so that 4iV le p tons is the actual 

ffhc 

number of leptonic degrees of freedom), and let N species = | J2 a (Qa/e) 2 denote the sum over 
all Dirac fermion fields weighted by the square of their electric charge assignments. Then, 



using the definition fl2.21|) of Q[x\ and the leading-log forms ( |3.28|) and ( |3.35| ) for the relevant 



contributions to the collision integral, the explicit form of the functional Q[x] becomes 

f^ 1 = f dp / »[i - m\ p 2 x e \p) 

JVleptons 7T Jo 

_ (e4lne _ 1} A|™ j™ dpm[ i- m] {p 2 [x e+ (p)T + 2x e+ (p) 2 } 

- (e 4 lne" 1 ) JL dppf^p) [1 + / »] X e \p? ■ (4.3) 

Varying the above leading- log approximation to Q[x e+ ] with respect to x e+ generates 
an ordinary differential equation for x e+ , as was noted in the context of shear viscosity 



23 Alternatively, 5j(p) is CP even, so Xi(p) = PiX(p) is CP even, which means xip) itself is CP odd. 
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by Heiselberg p8| . One option would be to solve that differential equation numerically. 
However, we find it both numerically and conceptually more convenient to instead directly 
attack the variational problem itself. This is also good practice for going beyond leading-log 
order, where the variational problem does not reduce to simple differential equations. 

Since every term in this leading log form for Q[x e+ ] is a one dimensional integral over 
p, it is possible to perform the maximization directly by discretizing values of p and then 
maximizing the resulting discrete quadratic form depending on a finite set of values of 
X e+ {p) ■ This is completely straightforward, but requires a very fine discretization in order to 
determine Q max with high accuracy. Another approach, which is both efficient and remains 
practical when applied to the full leading order in g calculation, is to maximize Q[x e+ ] within 
a variational subspace given by the vector space spanned by a suitable set of basis functions, 
{0( m )(p)}. In other words, one uses an ansatz 

TV 

X e \p)= (4-4) 

m=l 

where N is the size of the basis set considered, and the coefficients variational 
parameters which will be tuned to maximize Q[x e+ ]- Inserting the ansatz ( |4.4| ) into the 
functional fl4.3|) produces an N- dimensional quadratic form, 

N N 

Q\{_Q"rri\\ ^ t ®"m S m ^ ^ t dm C mn (l n , (4-5) 
m=l m,n=l 

where the basis set components of the source vector, S m = (^Si,<p\ m ^ and of the linearized 
collision operator, C mn = , C<f)[ n ^ may be read off from the previous expression ([4.31). 
Maximizing Q is now a trivial linear algebra exercise which gives a = C~ l S and 



a 



2 Q^ x = U-S = lS T C-'S, (4.6) 



3 ^ max 3 - ~ 3 



where a = \\a m \\ and S = \\S m \\ are the iV-component coefficient and source vectors, re- 
spectively, in the chosen basis, and C = \\C mn \\ is the (truncated) collision matrix. Given 
a particular choice of basis functions, the individual components {S m } and {C mn } may be 
computed by numerical quadrature, and then a single N x N matrix inverse yields the 
conductivity via the result ([4.6|). 

It remains to choose a good family of basis functions for the variational ansatz. There is 
a surprisingly good single function "basis set," namely 4>^\p) = p/T. Its use permits one 
to find analytic expressions for Si and Cn, 



^l-parameter = ^ ^ ? 



9C(3) 



7T 



^-parameter = Uptons + £ ) In e" 1 ) . (4.8) 



When substituted into Eq. ([4.6|), this yields our previously quoted approximate result 



Eq. ( |1 . 1 1|) . This one parameter variational ansatz is surprisingly accurate for several reasons. 



First, if one uses Boltzmann statistics instead of the correct Bose or Fermi distributions, 
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FIG. 4. Comparison of the exact (leading- log) form for x e {p) to the one and six parameter ansatz results, 
for the case of an e + e~ plasma (left) and a three lepton and five quark plasma (right). The six parameter 
ansatz curve is essentially indistinguishable from the exact curve except for p < 0.1T. The dashed curves 
look exactly the same for both of our two ansatze. 

[so that all appearances of /o(l±/o) are replaced by e _P//T ], then this one-parameter ansatz 
turns out to be exact. For large momenta, p > T, this modification produces negligible 
change in the integrands appearing in (|4.3|) , but it seriously mangles the integrands at small 
momenta, p < T. However, the dominant contribution to the integrals in Eq. ( f4.3|) comes 
from p/T ~ 4, where the effects of quantum statistics are already rather small. Additionally, 
at smaller momenta the Boltzmann approximation turns out to over-estimate the contribu- 
tion of diagram (C) to the scattering integral, but to under-estimate the contributions of 
diagrams (D) and (E). Hence there is some cancellation when the contributions are com- 
parable. Finally, as in any variational approach, the error in the resulting conductivity 
scales as the square of the error in the function x e+ (p)- For these reasons, a well chosen one 
parameter ansatz does surprisingly well. 

We have investigated two larger basis sets, each of which contains x = p/T as a special 
case.0 One choice for an iV-element basis set is 



4> (m \p) 



(p/T)' 



N-l ' 



m 



N, 



(4.9) 



The second is 



p 

T ' 



0M(p) = £ e -p/M 



m 



2,...,N, 



(4.10) 



with 



24 The basis sets 
order calculations 



9h and (4.10) were actually selected in order to yield good results in our full leading- 
For that application, it is helpful to have a non-orthogonal basis of strictly positive 



functions. For a leading-log analysis, many other simple choices of basis set would be equally good. 
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(c 2 ,c 3 ,c 4 , •••) = (log-, log-, log-, log-, log-, log-, log-, •••)• (4.11) 

For each basis set, we find that using the first four basis functions determines the (leading- 
log) value for o with a relative error less than 10~ 5 , while six basis functions are good to 
better than one part in 10 6 . The determination of x e+ (p) itself is also excellent; we compare 
the "exact" function x e+ {p) [obtained from a discretization with very large N] to the 1 and 6 
parameter variational ansatze in Fig. |j. The error barely visible in the six parameter ansatz 
for x e+ at p ^ T/10 is quite irrelevant to the determination of the conductivity because this 
region of the integration domain contributes virtually nothing to the relevant integrals. One 
may show that the exact (leading-log) x e+ (p) function vanishes non-analytically as p — > 
but this behavior has not been built into our choice of basis sets. Our numerical results for 
the leading-log conductivity have already been presented in Table [I] of the summary. They 
are based on six test functions, and are accurate to the last digit shown. 

The best previous result in the literature for the leading-log electric conductivity in a 
relativistic plasma is that of Baym and Heiselberg ||20|| . Their analysis differs from ours 
in the following respects. First, they neglect diagrams (D) and (E), and the same-charge 
contributions (particle on particle or anti-particle on anti-particle) in diagram (C)q Second, 
they consider only the one parameter ansatz, x{p) = p/T, and they determine its magnitude 
by integrating both sides of the Boltzmann equation against p, rather than against px{p)- 
Their treatment is thus not a variational approximation, and consequently the result for the 
conductivity is linearly, not quadratically, sensitive to the error in the x{p)- If we drop the 
contribution of diagrams (D) and (E) from our analysis, our result for a is 7.6% lower than 
theirs; including all diagrams, for a pure e + e~ plasma, our result for a is just under half of 
theirs. 



V. FLAVOR DIFFUSION 

As discussed in section [HJ the determination of flavor diffusion constants exactly parallels 
the calculation of conductivity, except for the use of differing charge assignments appropriate 
for the globally conserved current of interest, and an overall factor of the corresponding 
charge susceptibility E a = dn a /dfi a . 

We will consider a theory with a simple gauge group, some set of active matter fields 
(fermions and/or scalars) in arbitrary representations of the gauge group, and no other 
significant interactions. Hence, the net number density of each matter species (or 'flavor') 
is a conserved charge which will relax diffusively. We will only consider flavors carried 
exclusively by fermions. The diffusion constant for the net number density of species a is, 
from Eq. ( |2.26| ), given by 



25 They argue that e~e~ — > e~e~ and e + e + — » e + e + processes do not change the net current and so do 
not affect the conductivity. The final inference is incorrect except in the one-parameter ansatz. Though 
these collisions do not change the net current, they do change the distribution of velocities, which can then 
indirectly affect the rate at which the current is changed by other processes (although numerically the effect 
is quite small). Specifically within the one-parameter ansatz, however, the e + and e~ distributions are 
individually just normal thermal equilibrium distributions boosted parallel and anti-parallel to the electric 
field; and so for this restricted ansatz these scattering processes have no effect at all. 
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D(l q Qmax , / 1 — a ; (^-l) 



= l,<?=<3a 



where the functional Q[x] is to be maximized in the I = 1 channel with charge ±1 assigned 
to all particles or antiparticles of species a, respectively, and to all other excitations. To 
be consistent, all particles and anti-particles of species a should then also be included in the 
number density n a used to define the susceptibility S a [from Eq. ( |2.28| )1, giving 

s -=k =ld «~ T2 - (5 - 2) 

Once again, charge conjugation (or CP) invariance implies that the particle and antiparticle 
departures from equilibrium are equal and opposite, x b — ~X b ^ while the gauge boson 
departure from equilibrium must vanish, x 9 — 0.0 

In the leading-log approximation, the gauge boson exchange contributions ( |3.28| ) to the 
collision operator are diagonal in particle species. [That is, there are no x a X b cross terms 
in ( |3.28| ) .] And when symmetry prevents the gauge bosons from having any departure from 
equilibrium, the leading- log annihilation/Compton scattering contributions (|3.35|) to the 
collision operator are also diagonal in particle species. Consequently, in the case at hand, 
where the driving term Si in the linearized Boltzmann equation is only non-zero for species a 
and x g vanishes, the only non-zero departure from equilibrium emerging from the linearized 
Boltzmann equation (or the maximization of Q[x\) will be for species a. 

More physically, the essential point is that in the theory under consideration, in leading- 
log approximation, a departure from equilibrium in the net number density for species a 
will relax diffusively without inducing a net number density in any other species. Thus, 
for the determination of the diffusion constant D a , the variational functional Q[x] m &y 
be expressed solely in terms of x° Using the leading-log forms (|3.28| ) and ( |3.35|) for the 



relevant contributions to the collision integral, and inserting the species a particle number 
as the relevant charge in the source term ( |2.23| ), the explicit form of the functional Q[x] 
becomes 

Q[x " ] ~ y r dpfS(p)[i-fS(p)}p 2 x a ( P ) 



d Ra n 2 Jo 

- (# 4 ln# 



r> rffh 

4, -1\ U Ra 



967T 3 

ii„ C RaP 



i, 

oo 



dp fSip) [1 - fS(p)) {P 2 [X a (p)f + 2 X a (p) 2 } 



(g*\ng- 1 ) ^ dppf£(p) [1 + f Q {p)] X a {pY ■ (5.3) 



Inserting a finite basis-set expansion for x a (p)i as in Eq. ( |4.4j ), extracting the resulting basis 
set components for the source vector {S m } and collision matrix {C mn }, and inverting the 
(truncated) collision matrix to determine Qmax, proceeds exactly as described in the previous 
section. 



26 This follows because the gauge boson distribution, which comes from the color-singlet part of its density 
matrix, is charge conjugation even, while the conserved charge of interest is charge conjugation odd. 
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Using the same one-term ansatz x a ^ p/T with the functional (|5.3|) leads to our pre- 
viously quoted analytic approximation ( |1.23|) . Specialized to an SU(3) gauge theory with 
Nf fundamental representation quarks, this reduces to Eq. ( |1.22| ). Maximizing the above 
Q[x a ] using six or more terms of either basis set (|4.9| ) or (|4.10| ), and the specific values 
C*f = 4/3, Tp = 1/2, and Ta = 3 appropriate for SU(3), leads to the numerical values shown 
in Table |ffl|. 

Four previous determinations of the QCD diffusion coefficient merit mention. The first is 



by Heiselberg [ 19| . Within the single parameter ansatz, x a (p) ^ p/T, his analysis coincides 
with ours, except that he leaves out diagrams (D) and (E). He makes an additional factor 
of 16/9 error in his final value for the diffusion constant, which is therefore just more than 
twice ours for the case Af = 6. We are unable to trace the origin of this error; equations 
(9)-(12) of his paper appear correct, but (13) is not. Joyce, Prokopec, and Turok JH]] also 
treat quark diffusion. In this paper they make a leading log treatment in which they only 
consider diagram (C), which means they only keep the NfTp part of the second term in 
Eq. (|5.3|) . Furthermore, they neglect the energy transfer of a collision, which effectively 
means that they drop the p 2 [x a {p)'] 2 term in the associated integral but keep the 2[x a (p)] 2 
term. Under this "approximation," it is not difficult to solve exactly for the resulting 
form for x{p) without any recourse to an ansatz. Their result for the diffusion constant is 
about 5% higher than our result would be if we only included diagram (C) [which means 
retaining only the middle term in the denominator of Eq. ( |1.22j )l. In fact diagram (C) is 
sub dominant, so their result is fairly far off; for a six quark plasma their result is 2.3 times 
ours, and for fewer quarks it is worse. They also treat lepton diffusion, making the same 
approximation and a corresponding sized error. Another treatment of the baryon number 
diffusion constant is given by Moore and Prokopec |22j| . Their analysis coincides with our 
treatment if we make the one parameter ansatz x a {p) p/T. They did include diagrams 
(B) through (E). However, they failed to include the (1/2) symmetry factor in diagram 
(D) and make an algebraic error in diagram (E), so their results for these diagrams are off 
by almost a factor of 2, and their result for the resulting diffusion constant is about 12% 
smaller than ours. The treatment in another paper of Joyce, Prokopec, and Turok [23| is 
similar to the treatment of Moore and Prokopec, but only includes diagram (C) and makes 
some errors in its evaluation. 



VI. SHEAR VISCOSITY 

As shown in section |H], the solution of the linearized Boltzmann equation needed to 
determine the shear viscosity is equivalent to the maximization of the functional 

Q\x\ = Xij) - \(xij, CXij) , (6.1) 

now specialized to i = 2, with S£(p) = -Tf£(p)[l ± f£(p)} |p| %(£), ^(p) = iij(p) X a (p), 
^ij(p) = yf (piPj ~ an d the linearized collision operator C defined by Eq. ( |2.24| ). The 

resulting viscosity is proportional to the maximal value of Q[x], 

V = Tt Qmax • (6.2) 



35 



+ 



The source term Sij is invariant under all flavor symmetries as well as charge conjuga- 
tion and CP. Hence, the departures from equilibrium which solve the linearized Boltzmann 
equation = C\ij (or equivalently, which maximize Q[x]) will also be invariant under 
these symmetries. So, for example, in QCD, every quark and anti-quark species will have 
the same departure from equilibrium which we will denote as x q \ an d all gluons will share 
another common departure \ 9 . But there is no reason for \ q to equal x 5 \ nor f° r either one 
to vanish. Hence, two distinct functions must be varied independently to find the maximum 

of Q[x)- 

Using the leading- log forms (3.28) and ( [3.35 ) for the relevant contributions to the colli- 
sion integral, and specializing to a QCD-like theory, the two terms in Q[x 9 ,X q ]y Eq. (|6.1| ), 
explicitly equal 

(Si... j ,Xi~ j ) =— J o dpp 3 {d A f 9 (p)[l+f 9 (p)}x 9 (p) + 2d F N i / 9 (p)[l-/o 9 (p)]x 9 (p)}, 

(6.3) 

and 

^gHng-^ = ^ (TA+iVf Tf) r dP fM [l+fm ^ [X9{P)T + 6 

(T A +iV f 7» <*p fl(p) [1-f^p)} [p 2 [x q ( P )f + 6 X q {p?) 

+ - r *P/o(P) [1+/o(p)1 [x*(p) - X 9 (p)f • (6-4) 

167T J JO L J 

Our approach for maximizing Q[x] is an obvious generalization of the previous treatment 
used for conductivity or diffusion. To carry out the maximization within a finite dimensional 
variational subspace, we expand each of the undetermined functions in a finite basis set, 

N N 

X 9 (p) =E« m <P {m) (p) , X q (p) = E a N +m <P (m \p) , (6.5) 

m=l m=l 

with coefficients {a m }, m — 1, . . . ,2N, that are independent variational parameters. The 
variational functions (p) we use for shear viscosity are p times those presented in Eq. ( fl.9| ) 
or (|4.10|) ) because the conserved charge in question is now proportional to p. Inserting 
Eq. ( |6.5|) into expressions Q6.3) and (|6.4|) , one may read off the basis-set components of the 
source vector S and truncated scattering matrix C, 

(^ij: X.ij) ^ ' a m S m , [Xiji ^Xij) ^ ] a m C mn d n . (6.6) 

m m,n 

As before, one has Q[{a m }] = a T S — \a x C a. At the maximum, the vector of coefficients 
a = C" 1 ^, and 

V = JE = TE a ■ § = 55 S T C- l S. (6.7) 

The only change from the previous cases is that a set of N basis functions {cp^} now 
generates a 2N dimensional linear algebra problem. Note that non-zero block off-diagonal 
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components of C mn (those with m<N<n or m>N>n) represent a cross-coupling 
between the quarks and gauge bosons. These matrix elements arise only from the last term 
in Eq. (|6.4|) , which was generated by pair annihilation and Compton scattering processes 
[diagrams (D) and (E)]. If those contributions were absent, then Q would split into separate 
boson and fermion pieces, and the viscosity would be a sum of two independent terms, one 
due to fermions and the other due to the gauge bosons. 

The natural one-function ansatz is now (^^(p) = p 2 /T. Once again, this exactly solves 
the Boltzmann equation if the /o[l ± /o] quantum statistics factors are replaced by classical 
Boltzmann statistics. With this simple ansatz, one may evaluate the integrals in ( |6.3|) and 
(|6.4|) analytically. One finds 



S 



120C(5)T 3 



and 



C 



irT 3 



9d A 

x (g 4 ^g~ 1 



(d A C A + iV fC Z F C F ) 



7T Z 



d A C A 




d A 



(6.8) 



\N { d F C F 



971-2 

+ —N { d F C^d A 

LZo 



1 -1 
-1 1 



(6.9) 



Except for the numerical prefactors (and g 4 lng~ 1 ), these are exactly expressions ( 1.17Q and 
( 1.16| ) given in the Introduction. When combined with the relation (|6.7|) , these one-term 
ansatz results yield the approximate form ( |1.15|) quoted earlier. 

Maximizing Q[x g ,X 9 ] f° r an SU(3) theory, using six or more terms of a basis set with 



the form of either ( f4.9| ) or ( |4.10| ), but with each basis function multiplied by one additional 
factor of p, leads to the numerical values shown in Table [TJ. 

It should be noted that the role played by the final term of ( |6.4|) is not to relax the 
traceless part of the stress tensor (or momentum flux) of the plasma. At leading log order, 
the traceless stress is the same before and after an annihilation (or Compton scattering) pro- 
cess, because the outgoing momenta of the gauge bosons approximately equal the incoming 
momenta of the fermions. This is why the sum of these contributions vanish if \ g — X 9 - 
What these processes do is to transfer momentum between the fermions and the bosons. 
Such a transfer is important, because the gauge bosons and fermions equilibrate at different 
rates. The gauge bosons, with their larger gauge group Casimir, equilibrate faster than the 
fermions. Thus, the presence of a channel which transfers a departure from equilibrium in 
the slowly relaxing fermions, into a departure in the more rapidly relaxing bosons, speeds 
the relaxation process. 

This is why our result for 7] is somewhat lower than that found by Baym et al. |HJ , who 
missed diagrams (D) and (E), but whose treatment otherwise coincides with ours when we 
use the one term ansatz. For the case of three color QCD, the difference is largest at iVf = 2, 
where it is a 3% effect. Omitting these diagrams is significantly less important for viscosity 
than it is for conductivity and diffusion. This is partly because scattering via gauge boson 
exchange becomes more efficient at higher £ due to the £(£+1) factor in Eq. ( |3.28| ), and 
partly because the processes of diagrams (D) and (E) only indirectly affect the relaxation 
in this case, whereas for conductivity or diffusion the annihilation and Compton scattering 
processes can directly relax the relevant flux. 
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Abelian gauge theories are a special exception. In a QED-like theory, diagrams (A) and 
(B) are absent (because C\ = 0), and as a result diagrams (D) and (E) are essential. With- 
out these diagrams, photon fluctuations are not damped by any logarithmically enhanced 
process. If one omits these processes, then the matrix C becomes singular and the resulting 
leading-log viscosity diverges. This also helps explain why the QED viscosity we find in 
Eq. ( |1.18|) is so large. 



For a more complicated theory with a product gauge group, such as QCD plus QED, or 
the full standard model at temperatures above Myy, it is necessary to treat the departure 
from equilibrium of each species with differing gauge couplings as a distinct variable; so 
there are many independent x' s - There are also a large number of diagrams. However, if 
one of the gauge couplings is much smaller than the others, so that corrections of order of 
the fourth power of the ratio of couplings are negligible, then one may substantially simplify 
the calculation. This criterion certainly holds for QCD plus QED, where a;| M <C at 2 , and 
it is a reasonable approximation for the standard model, where {g'/g w ) 4 < 0.1. In this 
situation, all degrees of freedom which couple to the stronger gauge group, and all gauge 
bosons, may be regarded as remaining in equilibrium. Hence only weakly coupled fermions — 
leptons for QCD plus QED, right handed leptons for the high temperature standard model 
- are out of equilibrium. This approximation is valid for the more strongly interacting 
particles because their relaxation rates are much larger. In QCD plus QED, for example, a 
quark's departure from equilibrium relaxes at a rate which is 0[(q; s /q;em) 2 ] faster than the 
corresponding rate for the electron. This approximation is also valid for the more weakly 
interacting gauge bosons, such as the photon, because of the mixed weak/strong scattering 
processes represented by diagrams (D), (E), and (J), when the fermion is a quark, one gauge 
line is a gluon, and the other is a photon. This process occurs at a rate of order e 2 g 2 T, which 
is fast compared to the 0(e 4 T) (up to logs) relaxation rate of leptons. Similarly, in the hot 
standard model, hyper-photons scatter at a rate of order g g 2 T which is 0(a s /a') larger 
than the large angle scattering rates for right-handed leptons. This is the approximation we 
have used to obtain the viscosity results discussed in Sec. | for QCD plus QED, and for the 
standard model at high temperature. 

The shear viscosity in gauge theories has previously been considered by Baym, Monien, 



Pethick, and Ravenhall [ I7| . Their treatment is quite close to ours, and correctly analyzes 
the scattering contributions of diagrams (A)-(C) within the single function ansatz x cx p 2 . 
They also explore an ansatz with x oc p a in an effort to test the quality of the p 2 ansatz. 
However, they miss diagrams (D) and (E) which, as noted above, means that their values for 



rj are slightly too large. Shear viscosity has also been considered by Heiselberg fll8[ . For the 
case of pure glue (JVf = 0), Heiselberg gives a complete treatment, solving the variational 
problem for \ without recourse to any ansatz. We agree with his result for the Nf = 
viscosity, except for a trivial sign error in his presentation of the difference with respect to 
the p 2 ansatz value. However his treatment of the nonzero Nf case contains some errors not 
present in [1711, apparently from an incorrect treatment of statistical factors in diagram (B). 
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APPENDIX A: GROUP FACTOR NOTATION 



For any representation R of a simple gauge group, let {t a } denote the representation 
matrices of the Lie algebra generators. The normalization T R of that representation is 
defined by 

tr t a t b = T R 5 ab . (Al) 
The quadratic Casimir C R is defined by 

J2(t a ) 2 = C R l. (A2) 

a 

The two are related by 

T R d G = d R C R , (A3) 

where d R is the dimension of the representation R and dc=d A is the dimension of the gauge 
group. For the fundamental representation of SU(N), 

N 2 -l 

C F = ^ r , T F = l d F = N, (A4) 

while for the adjoint representation, 

C A = N, T A = N, d A = N 2 -l. (A5) 

The above notation is most natural to non-Abelian theories. However, results in the text 
can be translated to Abelian theories by taking = 1. Then, for each species c of charged 
particles, take dR — 1 and identify the generator t with the charge assignment (e c / e) of that 
species. The group factors are then 

C Rc = T Rc = e 2 Je\ C A = T A = 0, d F = 1. (A6) 
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